Interaction-induced dephasing of Aharonov-Bohm oscillations 
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We study the effect of the efectron-efectron interaction on the ampfitude of mesoscopic Aharonov- 
Bohm oscillations in quasi-one-dimensional (Q1D) diffusive rings. We show that the dephasing 
length Z/0 B governing the damping factor exp(— 2nR/ L^ B ) of the oscillations is parametrically 
different from the common dephasing length for the Q1D geometry. This is due to the fact that the 
dephasing is governed by energy transfers determined by the ring circumference 27ri?, making 
independent. 

PACS numbers: 73.23.-b, 73.63.Nm, 73.20.Fz 
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The Aharonov-Bohm (AB) oscillations of conductance 
are one of the most remarkable manifestations of electron 
phase coherence in mesoscopic samples. Quantum inter- 
ference of contributions of different electron paths in a 
ring threaded by a magnetic flux <E> makes the conduc- 
tance g an oscillatory function of <£, with a period of the 
flux quantum $o = hc/e; see Refs.[l||3|3jfor reviews. In a 
diffusive ring these $o _ P er iodic conductance oscillations 
are sample-specific (and would vanish upon the ensemble 
averaging), due to a random phase associated with diffu- 
sive paths. In this respect, the $o-periodic AB effect is 
a close relative of mesoscopic conductance fluctuations. 

Another type of the AB effect is induced by interfer- 
ence of time-reversed paths encircling the ring and is inti- 
mately related to the weak localization (WL) correction—. 
Its principal period is $0/2. It survives the ensemble 
averaging but is suppressed by a magnetic field pene- 
trating the sample. Below we concentrate on the first 
(mesoscopic, or $ -periodic) AB effect. Our results are, 
however, applicable to the second (weak-localization, or 
4>o/2-periodic) type of AB oscillations as well, as we dis- 
cuss at the end of the paper. 

In the present paper we will investigate the effect of 
the electron-electron interaction on the amplitude of the 
AB oscillations. Interaction-induced inelastic processes 
lead to dephasing of electrons, and thus to a damping of 
interference phenomena, in particular, of AB oscillations. 
The mesoscopic AB oscillations can thus serve as a "mea- 
suring device" for the electron decoherence. This idea 
was, in particular, implemented in recent experiments^^ 
where the low-temperature behavior of the dephasing 
time T<j, was studied, and two mechanisms of decoher- 
ence were identified: scattering off magnetic impurities* 
and electron-electron scattering. 

Let g(&) be the dimensionless (measured in units of 
e /h) conductance of the ring. It is convenient to expand 
the conductance fluctuations 6g($>) as a Fourier series 
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Within the conventional approach, when the dephasing 
time 1/t0 is introduced as a mass of the diffuson and 



cooperon propagators, Vd,c(q> <*>) ~ ^/(Dq 2 — iuj + 

the variance of the nth harmonic of the oscillations is 

suppressed by the factor 2 
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where = (Dt^) ' is the dephasing length, = 
(D/T) 1 / 2 the thermal length, D the diffusion constant, 
T the temperature, R the radius of the ring, and we set 
H = 1. [It is assumed in Eq. (0) that L^,Lt <C 2-kR] 
For a thin ring, is then identified with the dephas- 
ing length governing the WL correction in the quasi-one- 
dimensional (Q1D) geometry, which was found by Alt- 
shuler, Aronov, and KhmelnitskiiL^ to be 
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In fact, Aleiner and Blanter showed recently that rel- 
evant to the mesoscopic conductance fluctuations in wires 
has indeed the same form, Eq. Q , as the WL dephasing 
time. This seems to support the assumption that the de- 
phasing times governing different mesoscopic phenomena 
are identical. Equations J5J and © are commonly used 
for the analysis of experimental data. 

We will show below, however, that contrary to the 
naive expectations the formulas J3J and do not de- 
scribe correctly the dephasing of AB oscillations. Specifi- 
cally, if the interaction-induced exponential damping fac- 
tor of AB oscillations is presented in the form (<5<? 2 ) ~ 
exp(— 2tt Rn / L^ B ) , the corresponding length L^ B is 

parametrically different from Eq. ®. Moreover, L^ B 
depends on the system size R. 

We consider a ring coupled symmetrically by two leads 
to the bulk electrodes. We will assume the ring and the 
leads to be of a Q1D geometry, i.e. the widths of the 
wires are much larger than the Fermi wavelength but 
much smaller than R and . The only geometric pa- 
rameter characterizing the problem is then the ratio 7 of 
the resistance of the ring itself to the total resistance of 
the ring with the l eads . By definition, < 7 < 1. 

Following Refs. I7l8l l9lll(i the electron-electron inter- 
action can be represented by external time-dependent 



2 



random fields ip a (r, t) , with the correlation function 
(f a (r, t) <pP(r', t')) determined from the fluctuation- 
dissipation theorem, 



<^(r)^(r')> 6 



-Im U(r, y'-u) 5 aP coth— . (4) 



The conventional form for the dynamically screened 
Coulomb interaction in a diffusive system is& 
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where Uo(q) is the bare Coulomb interaction, H(q,u) — 
vDq 2 1 (Dq 2 — iu>) is the polarization operator, and v is 
the density of states. For a diffusive system, the bare 
Coulomb interaction may be neglected compared to the 
polarization since the inverse screening length, which is 
of the order of the Fermi wave number, is much larger 
than typical momenta q. As we will see below, the char- 
acteristic momenta q for our problem are of the order of 
the inverse system size R . In view of the nontrivial 
geometry of our system, it is thus more appropriate to 
work in the coordinate representation. A corresponding 
generalization of Eq. JSJ can be readily obtained, yielding 



Im^r/iwl-Imn-^r.r'iw) = ^(r,r') , (6) 

where T> is the propagator for the Laplace equation, 
— V 2 £>(r, r') = 8(r — r') with zero boundary conditions 
at the contacts with bulk electrodes. Substituting Eq. © 
in Eq. Q}, we get, for relevant frequencies w«T, 

2T 

<^(r, t) </(r', 0) = ^ V(r, r') 8 a[) S(t - t') . (7) 

Since the ring we are considering consists of Qf D wires, 
the corresponding diffusion propagator satisfies the one- 
dimensional diffusion equation 

[d t - Dd 2 x + i [<p° (x, t) - <pf>(x, t)} } V™g(x, t; x', t') 

= 5{x-x')5{t-t') (8) 

supplemented by appropriate matching conditions at 
junctions of the ring and leads. Here <5$ = <i>i — $2 
is difference in the AB flux between the two measure- 
ments, which is incorporated in the matching conditions. 
The conductance correlation function is given by the con- 
ventional two-diffuson diagrams yielding (we drop 
the prefactor of order unity) 
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where angular brackets denote averaging over the exter- 
nal fields, 5(t) is given by 



S(t) = 12ttT J U/ , (e 1 )/'( £a )^-^ 
= 3 7rT 3 t 2 sinh -2 (nTt) , 



and /(e) is the Fermi distribution function. The function 
S(t) is peaked at t = with a width T _1 . We will replace 
it below by the delta function 5(t). This is justified if the 
dephasing effect during the time t—t' ~ T _1 is negligible, 
i.e. (0 Q (x,t)0 Q (a;,t))r- 1 < 1. Using Eq. ©, we find 
that the latter condition is equivalent to the requirement 
that the conductance of the sample is much larger than 
the conductance quantum e 2 /h ~ (25 kSl) -1 . This condi- 
tion is well satisfied in typical experiments with metallic 
rings, thus justifying the replacement S(t) — ► S(t). 

We now express the diffusion propagators in Eq. (jSJ as 
path integrals. We are interested in the regime of strong 
dephasing, when the relevant paths propagate only inside 
the ring and do not extend into the leads (see below), ft 
is convenient to introduce the angular coordinate 9 on 
the ring (— 7r < 9 < n), with 9 = ±ir/2 corresponding to 
the junctions with the leads (Fig.[TJ. Expanding the con- 




FIG. 1: The sample geometry and the angular coordinate 
8 introduced above. The paths representing the saddle-point 
solution are shown. The geometric parameter 7 is defined as 
the ratio of the resistance of the ring without the leads to the 
resistance of the total sample, 7 = Ri/(R\ + R2) . 



ductance fluctuations in Fourier harmonics with respect 
to the flux, 5g($) — ► 5g n , we then get 
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where the path integral goes over pairs of paths 
9\{t), 9i (t) which have a relative winding number n. The 
"potential" V(di,6 2 ) in Eq. (TTJ is given by V (61,82) = 
{[4> a (9i) — Q (#2)] 2 }; its explicit form can be straight- 
forwardly obtained according to Eq. Q by solving the 
diffusion equation in the ring with leadsii We will only 
need below the form of V(8i, #2) for both coordinates be- 
ing in the same arm of the ring. For \9i\ < 7r/2 we find 



x~5(t-t')(Vll(x u x 2 ,t)Vjl(x 2 ,x u t')), (9) V(6i,8 2 ) = 
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the expression for > tt/2 follows from symmetry con- 
siderations. 

We consider first the fundamental harmonic (n = 1) of 
the AB oscillations; a generalization to higher harmon- 
ics, n = 2,3,... will be done in the end. For n — 1 the 
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relevant pairs of paths interfere after half encircling the 
ring in the opposite directions. We are interested in the 
regime of a relatively high temperature, when the de- 
phasing effect is strong. In this case, the path integral in 
Eq. can be evaluated via the saddle-point method. 
As has been mentioned above, the paths representing the 
saddle-point solution (instanton) do not extend into the 
leads. Indeed, exploring a part of a lead and returning 
back into the ring would only increase the actionof the 
path. It is clear from the symmetry considerations that 
the optimal paths satisfy 6\{t) = —02(t). Furthermore, it 
is easy to see that the optimal initial and final points are 
located at maximum distance from the leads, i.e. ©i = 
and 02 = 7r (or vice versa). To within exponential accu- 
racy, the problem is then reduced to that of a particle of 
mass R 2 / D tunneling with zero energy in the potential 
(Fig.0 



ATR 



< 9 < § 

2 ] > f < & < 7T 

(13) 

from 9 = to 9 — tt. Since the potential is composed 
of quadratic parts, the corresponding instanton action is 
easily calculated, yielding {8g\) oc e~ s with 



s = a 



T 1 / 2 R^ 2 
v x l 2 D 
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Here C 7 is a coefficient of order unity depending on the 
geometrical factor 7, 
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[2(7+1). 
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27 (l - 7 2 ) +7r + 2arcsin7 



C 7 is equal to 7r 5 / 2 /2 3 / 2 in the limit of long leads (7 — > 0) 
and to 7r 5 / 2 /4 in the limit of short leads (7 — » 1). The 



To calculate the preexponential factor, we have to take 
into account small fluctuations of the initial and final 
points 0i, 82 around their optimal values, as well as fluc- 
tuations of the paths 9\ (t) , 9i (t) around the instanton so- 
lution. The corresponding calculation is sketched in the 
Appendix. The final result for the variance of the har- 
monics of the mesoscopic AB oscillations reads (up to a 
numerical prefactor) 
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where n = 1,2,..., and the action S is given by Eq. i|14fl . 

Let us discuss the obtained result jTSjl , ijTil) . First of 
all, it is essentially different from what one would obtain 
by using the formulas (J2J, (J3J. Indeed, the exponent in 
Eq. (|16fl scales in a different way with the temperature 
and with the system size, as compared to Eqs. |PJ and 
© . It is instructive to write Eq. l(Tfi|l in a form analogous 
toEq. @, 
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thus defining the Aharonov-Bohm dephasing length L 
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which is much shorter than in the strong-dephasing 
regime L^ B <C 2irR . The corresponding dephasing rate 
1/ T £ B = D/(L§ B ) 2 is thus given by 
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FIG. 2: The potential V(9) of the tunneling problem, 
Eq. I13L plotted in units of ATR/(vD) for different values 
of 7. 



above calculation can be straightforwardly generalized 

to higher harmonics of the AB oscillations, n = 2, 3, 

The optimal paths still begin at = or tt but now per- 
form n/2 windings in the opposite directions. Therefore, 
the corresponding action is S n = nS. 



To shed more light on the physical reason for the dif- 
ference between the conventional Q1D formula and 
our result l(T%|l . (|T$|) . the following qualitative argument 
is instructive. Calculating perturbatively the dephasing 
rate using the formula © for the screened Coulomb in- 
teraction in a diffusive system, one gets 



dq T 



2tt vDq 2 ' 



(20) 



In the calculation of the dephasing rate in a wir o 8 i 9 i 10 , 
the arising infrared divergence is cut off self-consistcntly, 
since only processes with energy transfers oj > t^ 1 con- 
tribute. As a result, the lower limit of integration in 
Eq. (1201 is — LT 1 , yielding the result JSJl. On the 
other hand, in the case of the Aharonov-Bohm dephas- 
ing rate, the relevant paths have to encircle the ring. 
Therefore, despite the fact that L^ B <§; 2ttR, the low- 
momentum cutoff in Eq. H2(J|) is set by the inverse system 
size (2T7R)- 1 . This yields l/r^ B ~ TR/{vD), reproduc- 
ing (up to a numerical coefficient) the result (|19fl . 
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It is worth emphasizing that our result (|19|) for the 
dephasing rate depends also on the geometry of the leads 
through the coefficient C 7 . We note a certain similarity 
to the dependence of the dephasing rate in a ballistic AB- 
ring on the probe configuration recently found in Ref. flol 

As has been mentioned in the introduction, our re- 
sults are also applicable to the WL (/i/2e-periodic) AB 
oscillations. Their nth harmonic Sg^ L is determined 
by cooperon paths with winding number n. Assuming 
that the magnetic flux penetrating the sample is negligi- 
ble and comparing the path-integral representations for 
{5g 2 n ) and 5g™ L we find 



3TL- 



■Sg. 



WL 
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where L = irR/j. This implies, in particular, that the 
dephasing length L^ B , Eq. (JTSJ, is the same for both 
types of the AB effect. Equation H21(l is a generalization 
of the relation^ between the WL correction and conduc- 
tance fluctuations for single-connected geometries. 

To summarize, we have studied how the Aharonov- 
Bohm oscillations are suppressed by dephasing caused by 
the electron-electron interaction. Using the path-integral 
formalism and the instanton method, we have obtained 
the result (|16|l . I|14|l which is parametrically different 
from the naive expectation ©. This demonstrates 
that the AB dephasing rate 1 / r^ B , Eq. 119|) , is paramet- 
rically different from the dephasing rate I/t^, Eq. J3J), 
corresponding to a single-connected geometry. Physi- 
cally, the difference can be traced back to the fact that 
1 /Tfj, is determined self-consistently by the processes with 
energy transfers of the order of itself (or equivalently 
with momentum transfers ~ 1/L^), while the character- 
istic energy and momentum transfers governing l/r^ B 
are determined by the system size. For this reason, the 
Aharonov-Bohm dephasing rate l/r^ B depends on the 
ring radius R, diverging in the limit R — > oo . 

Valuable discussions with I.L. Aleiner, B.L. Altshuler, 
N.O. Birge, H. Bouchiat, M. Buttiker, V.I. Falko, and 
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APPENDIX 

To find the preexponential factor in (| 1 <j|) for the vari- 
ance of the amplitude of the AB oscillations, we have 
to perform several Gaussian integrations over deviations 
from the instanton solution in the path integral . We 
will only calculate the parametric dependence of the pre- 
factor, neglecting numerical factors of order unity. 

First, let us consider small offsets of the initial and 
final points of the paths from their optimal position. 
The second-order variation of the action S 2 S will be a 
quadratic form of the offsets, 5 2 S — Uij SQiSQj, where 
2 = 1,2. Using that 8 2 S ~ 1 for SOi ~ 1, we get 

(dettty)~ 1/2 ~ S- 1 . (A.l) 
Second, we have to account for small deviations of the 
paths from the instanton solution. The corresponding 
factor can be identified as the propagator for a harmonic 
oscillator with the parameters m ~ R 2 / D and moj 2 ~ 
RT/Dv. There are two such factors (one for each of the 
paths), yielding together 
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Finally, there is a Gaussian integration over the devia- 
tions of the time t spent on the path from its optimal 
value i opt ~ {vR/T) 1 / 2 . The corresponding factor can 
be estimated as 



d 2 S 
dt 2 
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Combining Eqs. ifTT]). (|A.ljl . I|A.2|) and (|A.3jl . we obtain 
the total preexponential factor as given in Eq. (|16fl . 
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